Dimensionality effects in dipolar fluids 
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Using classical density functional theory (DFT) in a modified mean-field approximation we in- 
vestigate the fluid phase behavior of quasi-two dimensional dipolar fluids confined to a plane. The 
particles carry three-dimensional dipole moments and interact via a combination of hard-sphere, 
van-der-Waals, and dipolar interactions. The DFT predicts complex phase behavior involving first- 
and second-order isotropic-to-ferroelectric transitions, where the ferroelectric ordering is character- 
ized by global polarization within the plane. We compare this phase behavior, particularly the onset 
of ferroelectric ordering and the related tricritical points, with corresponding three-dimensional sys- 
tems, slab-like systems (with finite extension into the third direction), and true two-dimensional 
systems with two-dimensional dipole moments. 

PACS numbers: 68.03. Fg, 64.70.Ja, 64.75.Xc 



I. INTRODUCTION 

Two-dimensional (2D) fluids consisting of parti- 
cles with classical dipole-dipole interactions such as 
(para)magnetic nanoparticles at interfaces [2-01, cobalt 
nanocrystals on solid surfaces Q , and suspensions of po- 
larizable colloids in 2D dielectrophoretic set-ups (EH, 
currently attract much attention. Indeed, as a result of 
the directionality of the interactions whose details can 
be tuned by external (in-plane, out-of-plane, or tilted) 
fields, 2D dipolar systems display a variety of interesting 
structures such as chains and bundles at low densities 
[2, [3 but also various solid phases [jjj]. Especially the 
self-assembled low-density structures suggest that such 
systems are promising candidates as tunable advanced 
materials Q with applications in electrical engineering, 
sensors [|[ and molecular miniature devices. 

For theory and computer simulations, exploring the 
full structural and phase behavior of 2D dipolar systems 
remains challenging. Apart from the above-mentioned 
aggregation phenomena, one topic investigated particu- 
larly by computer simulations concerns the appearance 
and characteristics of vapor-liquid transitions |9l-fl4| . An- 
other question touches the structure at high densities 
close to the range where crystallization is expected to 
occur. Various Monte Carlo (MC) simulation studies 
[t| [l5[ revealed the appearance of ferroelectric (or ferro- 
magnetic, respectively) domains, but overall frustrated 
(vortex) structures without true long-range orientational 
ordering. This is consistent with integral equation results 

lEBli where predictions on the low-temperature behav- 
ior are extracted by analyzing correlation functions. On 
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the other hand, recent Molecular Dynamics (MD) simu- 
lations [l3| revealed long-range ferroelectric ordering in 
dense, 2D Stockmayer fluids, where the dipole-dipole in- 
teractions are supplemented by isotropic Lennard- Jones 
(LJ) interactions. 

Similar to the dense fluid state, the nature of the 
2D crystalline structures formed at finite temperatures 
remains so far unclear [3, although ground state 
calculations indicate ferromagnetism for certain 2D lat- 
tice types such as hexagonal lattices [2(|. In three- 
dimensional (3D) systems the existence of long-range fer- 
roelectric (-magnetic) ordering under appropriate bound- 
ary conditions is well established [21-23] . Moreover, com- 
puter simulations of slab-like systems [24| , where the par- 
ticles are confined between two plane-parallel walls, have 
indicated that this type of confinement can actually pro- 
mote long-range ordering of the dipole moments, if the 
wall separation L z is sufficiently large. However, decreas- 
ing L z to values where less than three monolayers can 
form, the ordering seems to disappear [25j j . 

The diversity of simulation results shows that spatial 
dimension has a profound influence of the ordering be- 
havior of dipolar fluids. The purpose of the present study 
is to collect and compare theoretical results on that is- 
sue based on a relatively simple, mean-field like approach. 
Specifically, we employ classical density functional th eory 
(DFT) in the modified mean field approximation [26- 28] , 
where the pair correlations are replaced by their low- 
density limit, i.e., the Boltzmann factor. The applica- 
tion of this approach for three-dimensional (3D) dipo- 
lar systems and their mixtures [29j was put forward by 
Groh and Dietrich [30l - l33l |. who considered Stockmayer 
fluids. Later the modified mean-field DFT approach has 
been used to study confined, slab-like Stockmayer fluids 
[34l |35| , with different degrees of sophistication regarding 
the hard-sphere part of the density functional, 

Here we apply the approach to a quasi-2D dipolar 
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(Stockmayer) fluid, where the particles are confined to 
a plane, but carry 3D dipole moments. Evaluating the 
phase diagram and comparing with corresponding DFT 
results for 3D systems, slab-like systems, and true 2D sys- 
tems with 2D dipole moments we can identify, on a mean- 
field level, the influence of spatial dimension and of the 
dimension of the order parameter on global ordering in 
fluid-like dipolar systems. Based on previous experiences 
one would expect that the mean-field approach for the 
quasi-2D system will (as it generally does) overestimate 
the stability of orientationally ordered phases. However, 
given the importance of meanficld-like approaches in the 
general context of spin and dipolar systems, and realizing 
that the mean-field DFT approach is, so far, still the only 
theory targeting the whole (homogeneous) phase diagram 
of dipolar systems, we think our results are important for 
a complete understanding of such systems. 

The remainder of the paper is organized as follows. In 
Sec. [II] we formulate the quasi-2D model and briefly de- 
tail the derivation of the corresponding grand-canonical 
functional. Numerical results for the phase diagram at a 
typical dipole moment are presented in Sec. IIIII There 
we also use Landau expansions to compare the onset of 
ordering in the quasi-2D system with the cases of 3D, 
slab-like and true 2D systems. Finally, in Sec. IIVI we 
summarize our results. 



II. THEORY 

The quasi-2D Stockmayer fluid consists of disk-like 
particles of diameter or at positions = (xj,i/j) in the 
x-y plane. The orientation of their 3D dipole moments 
fit is represented by the Euler angles = (6>i, </>;). The 
microscopic interactions between the particles stem from 
anisotropic dipole-dipole and isotropic LJ forces. The 
resulting pair potential between two particles with coor- 
dinates (1) = (ri,tJi) and (2) = (r 2 ,W2) is given as 



u(l,2) = 



00, r i2 < cr T 

Udi P (ri2, u 2 ) + u hJ {r\ 2 ), ri2 > cr T , 



(1) 



where ri 2 = ri — r 2 is the connecting vector be- 
tween the two particles, and r\ 2 — \c\ 2 \. Further, 



«dip(ri2,wi,<*>2) = {p 2 /rl 2 
? 12) (£2(^2) • £12) 



£i(wi) • £2(^2) ~ 3(£i(wi) 



is the 3D dipole-dipole interaction 

potential, and WLi{r 12 ) = 4e {a/r 12 ^j - (a/rxz) is 
the LJ potential. To mimic the fact that the effec- 
tive range of the r~ 12 repulsion varies with the ther- 
modynamical parameters, we choose in our DFT cal- 
culations a temperature-dependent hard core defined 
via the Barker- Henderson formula [36J, that is, ctt = 
f° dr (1 — exp [— /3itLj(r)]), where ft — l/k^T (with fee 
and T being the Boltzmann constant and the tempera- 
ture, respectively). 

To analyze the phase behavior we employ classical 
DFT, where the key quantity is the grand canonical po- 



tential fl as a functional of the singlet density p(r, oS) — 
^ SiLi <5( r— Yi)5(uj— Wj)^, with N being the total number 
of particles (37} • We restrict the analysis to fluid- like but 
possibly polarized ordered phases of the quasi-2D Stock- 
mayer fluid (in the following we assume, without loss of 
generality, the dipoles to be of electric nature) . In prin- 
ciple, investigation of this situation requires to perform 
a free minimization for the profile p(r,u)), thereby al- 
lowing the system to form domains (or other patterns 
with spatially varying polarization). Practically, how- 
ever, minimization including pattern formation is a quite 
challenging task as demonstrated in [32|, [33[. In the 
present study, where we are interested in the general ten- 
dency for ordering, we neglect that problem. That is, 
we focus on the polarization within a (macroscopically 
large) fluid domain, which may be part of a globally 
unpolarized system. We thus consider the singlet den- 
sity p(r,u>) — pa{uS), with p being the constant number 
density of particles and a(u>) the orientational distribu- 
tion function of their dipole moments. This function is 
normalized to I (i.e^, J a.(uS)duj — 1) and equals 1/4tt 
for isotropic states 34] . To describe an orientational or- 
dering along a specific direction, we expand the orienta- 
tional distribution function in terms of spherical harmon- 
ics, that is, a(ui) = EniEnp-i^m^raH] witri the 
coefficients a; m representing orientational order param- 
eters pjij . Nonzero order parameters with I = 1 corre- 
spond to a macroscopic polarization P = pp,p, with p = 
J a(ui)fi.(u))du). Indeed, transforming the Cartesian com- 
ponents of p = (Px,Py,Pz) in terms of spherical harmon- 
ics one obtains p = (— 2Re(an), 2Im(o!ix), ccio) 
[3~i] | . Similarly, nonzero order parameters ai m with I = 2 
indicate that there is a preferred orientation of the dipole 
axes. 



Within the grand canonical formalism, the system is 
characterized by its size A, the chemical potential /x c hem 
and the inverse temperature ft. In the present study, 
the dipolar contribution to the excess (interaction) free 
energy is treated in the modified mean-field approxima- 
tion, where the pair correlations are approximated by the 
Boltzmann factor 428j . In addition, following our pre- 
vious study on slab- like systems [34| . we perform a per- 
turbation expansion of the Mayer function and truncate 
the latter after the quadratic term. Such a truncation has 
been first employed in a DFT study of the surface tension 
of polar fluids [261 ] . Later, results for the phase diagrams 
of bulk polar fluids [3l| have indicated that the second- 
order theory yields data very close to those from the full 
modified meanfield approximation (without any trunca- 
tion) . As a consequence of the truncation, the resulting 
excess free energy contains only terms up to I = 2. For 
a detailed calculation for the quasi- 2D case (which pro- 
ceeds analogous to the slab case) we refer the reader to 
Ref . [38j . The resulting expression for the grand canoni- 
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cal functional is given by 



a) 



^ = P fin (pA 2 ) - 1 - /3/i c hem + J dua(uj) In [47ra(o;)]^ 



pln(l - rjr) + p 

7rp 2 a 2 9l (T) 
4ir 2 



1 7]T 



+ —p 2 p 2 ^ 1 +n~ 1 g 2 {T)){a 2 w ~\a 11 \ 2 ) 
H A p 2 l3 2 a-%(T)G[a(uj)}. 



On the right side of Eq. © , the first line contains ideal 
gas contributions (involving the thermal wavelength A = 
h/ y/2ir/3m) and the orientational entropy (last term). 
The second line contains the excess free energy of our 
reference system, the hard disk fluid 39], involving the 
2D packing fraction t\t = Trpa^/A. The three last terms 
account for the dipolar interactions, where the func- 



tions gi(T) = 
9*(X) = /„7t) 



= X^t) d x x (exp [— 4/3e(x 
dxx~ 2 (exp \—Aj3e(x~ 12 - 



12 



')] 



3 )] 
1) 



and 



9s(T) = J™ T) dxx~ 5 (exp[-4l3e(x~ 12 - x~ 6 )]), with 
a(T) — (Tx/cr. Finally, in the last term on the right 
side of Eq. ©, G[a(uj)] = ha 2 + t 2 a oa 2 o 
here t\ = 
3tt 2 /(15V2), U 



_ halo + 

U\a 21 \ 2 + t 5 \a 2 2\ 2 , where tj = 10tt 2 /(15V2), t 2 = 



-2^/(15^2) 



-2V5tt 2 /(15V2), ts 
and t 5 = n 2 /(15V2). 

Minimization of the functional ([2]) with respect to p 
and a(uj) yields the Euler-Lagrange equations for this 
problem [38|]. They consist of a set of nonlinear, coupled 
equations for the density p and the OPs ai m appearing 
in Eq. ([2]). The equations are solved numerically using a 
Newton- Raphson algorithm [40| . 

In the following we characterize the state of the 
quasi-2D Stockmayer fluid by the dimensionless density 
p* = per 2 , temperature T* — k^T/c, chemical potential 
Mchcm = e_1 (Mchcm - 2fc B Tln( A/cr) ) , and dipole moment 
m* = p,/ \J ea 3 . The quantity to* 2 measures the strength 
of the dipolar interactions in an antiparallel side-by-side 
configuration relative to the LJ interactions. We note 
that the coupling parameters to* and T* are equivalently 
defined in 3D (or slit-pore-) dipolar systems [3J], so that 
the results can be conveniently compared. 



III. RESULTS 

Following earlier DFT studies on confined Stockmayer 
fluids [34| we consider a system characterized by to* = 
1.5, a typical value for moderately polar molecular fluids 
such as chloroform [4l|. The calculated fluid phase dia- 
gram in the density-temperature and chemical potential- 
temperature plane is shown in Figs.[T^) andQ]}), respec- 
tively. The latter representation better relates to typical 
sorption experiments (42[. For small and intermediate 
densities (or chemical potentials) we find a state where 



(2) 



b) 




FIG. 1: (Color online) a) Phase diagram of the quasi-2D 
Stockmayer fluid (m* = 1.5) in the density-temperature 
plane. The regime under the black stars indicates the IF- 
FF coexistence region. The blue dashed line is the line of 
critical points [see Eq. Q], which starts at a tricritical point 
(TCP). Included are results (from 34]) for the TCPs of the 
corresponding 3D Stockmayer fluid (red diamonds), and two 
slit-pore systems with wall separations L z = 10a (black tri- 
angle) and L z = 4u (black circle). The red squares indicate 
the IG-IL coexistence curve of the quasi-2D system resulting 
from a separat calculation where the system is forced to be 
unpolarized. b) Corresponding phase diagram in the chemical 
potential-temperature plane. 



all order parameters ai m are equal to zero, and those 
with I = 2 are either zero (to = ±1,±2) or negative 
(to = 0). Thus, there is no global polarization and nei- 
ther a global ordering of the dipole axes; we therefore 
refer to this state as "isotropic fluid" (IF). The negative 
values of a 2 Q merely indicate that the dipoles tend to 
avoid to be oriented parallel or antiparallel to the z-axis; 
rather they prefer to lie (with random orientations) in the 
x-y-plane. This is an expected effect in a dilute, quasi- 
2D dipolar system (consistent with simulations and other 
theoretical studies, see e.g. [II]). At higher densities, the 
system then develops a non-zero polarization, to which 
we refer to as "ferroelectric fluid" (FF). Within this state, 
the vector P points along an (arbitrary) direction in the 
x-y-plane (as reflected by ai,±i ^ 0, aio = 0). Notice 
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that the preference of in-plane polarization (rather than 
out-of-plane polarization) can already be seen from the 
prefactor of the corresponding terms oc |ai m | 2 in Eq. ((2|). 

The transition between the IF and the FF phase is 
discontinuous in a; TO and p (yet not in ^ c hem) for tem- 
peratures below a tricritical temperature T t * = 1.57 [see 
Fig- UK)]- Above the tricritical point (TCP) the transi- 
tion becomes continuous, which results in a line of critical 
points. The appearance of a TCP is a typical feature of 
DFT predictions of the phase diagrams of dipolar flu- 
ids Ho, HH HI| and also Heisenberg fluids [43| . Recent 
MC studies for 3D dipolar fluids confirm that the transi- 
tion between isotropic and ferroelectric fluid is of second- 
order in a broad range of temperatures (23|. Within 
the DFT, the line of critical points can be determined 
from a Landau expansion of the free energy, assuming 
that the OPs characterizing the FF state are small (i.e., 
\ot\m\ <C aoo = l/\/47r). To this end we expand the in- 
tegral /[a(cj)] = J duja{<jj) In [47ra(w)] in Eq. @, that is, 
the orientational entropy, in a Taylor serios around the 
isotropic state (where a{uj) — a o^oo( w ) = 1/4tt). Col- 
lecting those terms in the resulting approximate func- 
tional, ft, that are proportional to ai ±1, we obtain (38j 

^ oc laufi^p-^p 2 ^ 1 + *- 1 92(T))), (3) 

where the first term stems from the orientational entropy, 
whereas the second term results from the interaction free 
energy in Eq. ([2]). The second order phase transition is 
characterized by a change of sign of the factor of |an| 2 
in Eq. ©. We thus obtain 

fc B T c — 2D = ^P^T 1 + v~ 1 92{T c )). (4) 

This is approximatively the equation of a straight line in 
the density-temperature plane, consistent with what one 
sees in Fig. QJ,). 

In Figs. Hk) and b) we have included DFT data for 
tricritical points of Stockmayer fluids in 3D and in slit- 
pore geometries. Within the latter situation, the par- 
ticles are confined between two planar, attractive walls 
separated by a distance L z [34|. We note that, both for 
the 3D and the slab case, our data somewhat differ nu- 
merically from those presented in another recent DFT 
study of confined Stockmayer fluids (35|. This is since 
we used (contrary to [Hj]) a temperature-dependent par- 
ticle diameter and a homogeneous ansatz for the num- 
ber density in the slit-pore. However, from a qualitative 
point of view, the observed trends regarding the impact 
of confinement on the TCP are the same in both studies 
[3~il [351 1 . In particular, both predict that decreasing L z 
(that is, increasing the degree of confinement) shifts the 
TCP towards lower temperatures and somewhat lower 
densities For example, according to [34], the tricritical 
parameters (p* cp ,T^ ) are (0.62,1.92) at L z — 00, i.e., 
in the 3D (bulk) limit, (0.60,1.80) at L z = 10a and 
(0.58,1.68) at L z = 4a 34]. Consistent with this ten- 
dency, the TCP of the quasi-2D Stockmayer fluid (which 



corresponds to the limit L z — > 0) is found at even lower 
density and temperature, specifically at p* cp — 0.52 and 
T* cp — 1.57. A somewhat different behavior emerges 
in the chemical potential-temperature representation de- 
picted in Fig. [TJd) . From the location of the TCPs, we 
see that the tricritical chemical potential, H^ em , of the 
quasi-2D fluid (L z — 0) is only slightly smaller than that 
of its 3D counterpart (L z — 00). On the other hand, the 
corresponding values for At c ^g m of confined Stockmayer 
fluids (L z = 10a, 4a) are significantly smaller (consistent 
with [35J). We attribute this non-monotonic behavior of 
Mchem u P on lowering of L z [see Fig.QJ))] to the fact that, 
for the confined Stockmayer fluids, the walls were con- 
sidered to be attractive [34]. This feature is known to 
support capillary condensation (or, more generally, the 
formation of denser phases), accompanied by lowering of 
the critical value of ^ c hcm- The particles in the quasi-2D 
system do not interact with any walls; thus, there is no 
capillary condensation phenomenon. As a result, Mchem 
for the quasi-2D system nearly agrees with the 3D value. 

The reduction of spatial dimension not only shifts the 
TCP, it also has a profound influence on the topology of 
the phase diagram. Indeed, while the 3D Stockmayer 
fluid with the same dipole moment (m* = 1.5) exhibits, 
in addition to the IF-FF transition, a condensation tran- 
sition within the isotropic liquid (IL) phase, such a tran- 
sition is absent in the quasi-2D system (see Fig. [TJ. We 
can artificially stabilize a condensation transition in the 
2D system by setting all order parameters (except from 
CK20) t° zero. The result of this calculation is indicated 
in Fig. [Th.) by the red squares. It turns out that the IG- 
IL critical point (T* = 0.98) is located within the IF-FF 
phase coexistence region; therefore it is thermodynami- 
cally unstable. However, such IL configurations may still 
be relevant in the context of the phase separation kinetics 
(i.e., in non-equilibrium situations), where they can oc- 
cur as transient states during the change from the IF to 
the FF state at a temperature T < T*. We also note that 
the suppression of the IG-IL critical point is consistent 
with previous DFT results (at m* = 1.5) in very narrow 
slit-pores, such as L z — 4a [34J. On the other hand, MC 
simulations for quasi-2D Stockmayer fluids predict stable 
isotropic liquid phases for dipole moments up to (at least) 
to* = y/6 [l3,0,[l7|]. Therefore, the DFT seems to over- 
estimate the stability of the ferroelectric phase, similar 
as it does in 3D [34| . Within the DFT, one would ex- 
pect a recovery of the isotropic liquid state when the LJ 
attraction finally dominates then dipolar coupling, i.e., 
when to* decreases towards even smaller values. 

We now discuss in more detail the influence of spa- 
tial dimension and its interplay with the dipolar cou- 
pling strength on the tricritical point (/9* cp , T t * cp ), above 
of which the low-temperature discontinuous transition 
between the IF and FF states changes into a (line of) 
second-order transition(s). Specifically, we are interested 
in the position of the TCPs, in the quasi-2D system and 
its 3D counterpart, as functions of the parameter to*. In 
previous DFT studies [3(J HH, HiJ it was already shown 
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that the coupling strength influences the quantity T* cp 
much more than jO* cp , at least as long as m* > 1.0. 
Therefore, to estimate the dependence of T t * cp on m* in 
the quasi-2D system, we set the density equal to tricrit- 
ical density at m* = 1.5, that is, to p* \ 

quasi-2D 

= 0.52. 

The resulting function T* cp (to*) can then be easily de- 
termined from Eq. The same procedure is used for 
the 3D case, where the analog of Eq. © reads 44 1 



••L> 167r 2 

oc |an| 2 (4^ b - — — p^pl) (5) 



V 

with V being the volume and ph = p*a 3 being the den- 
sity of the bulk system. Equation ([5]) yields kBT^ = 
(4ir/9)pbp 2 ■ Fixing the density to that of the tricritical 
point at to* = 1.5 [34], Pt cp \3D = 0.61, we can again 
estimate the function T t * (to*) for a range of dipole mo- 
ments. Numerical results for the quasi-2D system and its 
3D analog are plotted in Fig. [5] For both systems, the 
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FIG. 2: (Color online) Tricritical temperature T t * p as func- 
tion of the dipole moment m* in the quasi-2D, 3D, and true 
2D geometry. In all cases, the density has been set to that 
obtained at m* = 1.5. 

tricritical temperature increases with to*, as one may 
expect when the dipolar interactions (which stabilize the 
FF state) become more and more important as compared 
to the spherical attractive ones. More interestingly, Fig. [5] 
reveals that reduction of spatial dimension shifts the tri- 
critical temperatures towards lower values at fixed to*. 
This shift can be reasoned from Eqs. ([3]) and ([5]): in both 
the quasi-2D and the 3D system, ordering competes with 
the same amount of (orientational) entropy, but the as- 
sociated decrease of interaction energy (oc p 2 p 2 ) is less 
pronounced in the quasi-2D system. From a physical 
point of view, this diminishment is a consequence of the 
reduction of the number of neighbors in a 2D system rel- 
ative to the bulk case. We also note another interesting 
point: whereas in the 3D system, only the long-range 
dipolar interactions contribute to the onset of ordering 
[see Eq. §5$], the corresponding onset in the quasi-2D 
system is also affected by the short-range interactions, 



as reflected by the appearance of the function gi (T) in 
Eq. ([3]). Since the gi function increases upon cooling 
down, the curves T t * p (m*) in Fig. [5] for the quasi-2D 
and the 3D system, respectively, cross each other near 
to* = 1. However, since we assumed a constant tricriti- 
cal density, the precise position of this crossing in Fig. [5] 
should be considered with some caution. 

Finally, we briefly discuss the influence of the dimen- 
sion of the dipole vector (rather than that of the space 
accessible for the particles) on the appearance of fer- 
roelectric order. Specifically, we consider a "true" 2D 
system where, in addition to the spatial confinement of 
the particles within the x-y-plane, the orientations of the 
dipole vectors are restricted to that plane as well. Indeed, 
as shown in previous simulations and theoretical studies 
(see, e.g., [15|) of 2D systems with freely rotating (i.e., 
three-dimensional) dipoles, these have a strong tendency 
to tilt into the confining plane especially at large cou- 
pling strengths. The "true" 2D system is therefore not 
completely unphysical. In the true 2D case, the orien- 
tational distribution function a{ui) depends only on one 
angle, 0, which describes the orientation of p, t relative 
to, say, the £-axis. To obtain the grand canonical func- 
tional we expand a(<p) in a basis of exponential functions 
[lH, i.e., a(4>) oc J2 m a m exp [imcj)]. Performing then the 
same Landau expansion as for the quasi-2D system, and 
isolating the terms proportional to the polarization, i.e. 
to a m —±i, the analog of Eq. ((3j) reads 



■2D 



A 



cx \a 1 \ 2 (27Tp-7r 2 p 2 p 2 /3(a^ 1 +a- 1 g 2 (T))). (6) 



A direct comparison of Eq. ([6]) with its quasi-2D ana- 
log in Eq. ([3]) shows that, at fixed density, the ferroelec- 
tric ordering in the true 2D system occurs at a higher 
temperature. This is a consequence of the decrease of 
the dipolar fluctuations (and thus, the orientational en- 
tropy) due to their restriction to the plane. Moreover, as 
revealed in Fig. [2] by the corresponding curve T t * cp (to*), 
the ordering is even promoted relative to the 3D case. 
This is consistent with tendencies found in a recent in- 
tegral equation study [3], although the latter predicts, 
for the low-temperature behavior, large ferroelectric do- 
mains rather than true long-range ferroelectric order. 



IV. SUMMARY AND CONCLUSIONS 

In this work we have calculated the fluid phase dia- 
gram of a quasi-2D Stockmayer fluid by means of density 
functional theory in the modified mean-field approxima- 
tion. At the dipole moment considered (to* = 1.5) the 
system exhibits an isotropic fluid phase where the dipole 
moments are randomly oriented, yet with a preference for 
in-plane directions, and a ferroelectric fluid phase char- 
acterized by global, in-plane polarization. Apart from 
exploring the quasi-2D phase behavior, another focus of 
our study was to identify the role of the dimension of 
accessible space, as well as that of the dimension of the 
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dipole vector. To quantify these effects on a mean- field 
level, we have considered the location of the tricritical 
point. Regarding the impact of space dimension, we 
have found that decreasing the system's dimension in in- 
direction from the bulk limit (L z — > oo) over slab systems 
(L z = 10cr, 4cr) towards the 2D limit (L z — > 0) shifts the 
TCP towards lower temperatures and densities. Further- 
more, the disappearance of the isotropic liquid phase in 
the quasi-2D system also shows that the confinement en- 
hances the stability of dense ordered phases relative to 
disordered ones. Clearly, care has to be taken with re- 
spect to the predictions of our mcan-ficld-like DFT ap- 
proach on a quantitative level. Indeed, from computer 
simulations EH E3] it is known that a quasi-2D Stock- 
mayer fluid at m* =1.5 does have a stable isotropic liquid 
phase at densities beyond the isotropic vapor-liquid criti- 
cal point, which is absent in our study. This discrepancy 
reflects the well-known tendency of the DFT to overesti- 
mate the stability of ordered phases. However, based on 
previous DFT studies for bulk and confined systems one 
would expect a recovery of the isotropic liquid state in the 
quasi-2D case upon further decrease of in* . A further in- 
teresting result of our study concerns the role of the spin 
dimension. Here we have found that complete restriction 



of the dipole moments on in-plane directions yields fer- 
roelectric ordering at temperatures not only higher than 
those in the quasi-2D system, but even higher than those 
in 3D. 

There remains the question whether fluid states with 
long-range ferroelectric order, as predicted by DFT, exist 
at all in quasi-2D and true 2D systems. As mentioned 
in the introduction, computer simulations give conflict- 
ing answers, which may also depend on the number of 
particles considered in the simulation (indeed, the MD 
study on quasi-2D systems by Ouyang et al. [17j . which 
does predict long-range ferroelectric ordering, involves a 
rather small system size). We should therefore interpret 
the present DFT results, which rely on the assumption 
of a spatially homogeneous orientational structure, such 
that the 2D geometry definitely promotes the existence 
of large ferroelectric domains, but not necessarily true 
long-range order. 

Despite these pitfalls, the DFT approach provides a 
general overview of the phase diagrams and highlights 
the dimensionality effects by providing the leading order 
terms in the free energy. From that perspective, it would 
be interesting to extend the study towards 2D systems 
in external fields. 
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